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Abstract— Deep space exploration missions are aimed at acquiring information about the solar 


system and a significant communication capacity has to be planned to transfer data for such very 


large distances [1]. Terrestrial atmospheric impairments on the space-to-Earth propagating signal 
are the major responsible for the signal degradation thus reducing the channel temporal availability 


[2]-[4]. In this work weather forecast models, coupled with microphysically-oriented radio- 
propagation models, are described in order to evaluate atmospheric effects at Ka-band. Estimation 


data return techniques are summarized and numerical results in a simulated operational scenario 
are illustrated in terms of received data volume using the BepiColombo mission as a baseline 
example. 
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e 1. Introduction 


The exposed topics are developed in the framework of the 
RadioMeteorological Operations Planner (RadioMetOP) 
program taking The BepiColombo (BC) mission (the European 
Space Agency (ESA) cornerstone mission to Mercury that is 
expected to begin in 2016) as a baseline example. The ground 
segment of the BC mission will be the Cebreros receiving 
ground station in Spain. 


GOALS: 


* Enhancing the predictability of the Ka-band downlink 
operations 


e Increasing the received data volume at G/S 


* Exploit weather forecast to decide how to set up the 
downlink between spacecraft and Eart 





2. RadioMet Operations 
UE Planner blocks 
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° Xm 4D (space-time) meteorological state vector = 
thermodynamic variables + microphysical variables 

e Li, = path attenuation 

* Ti+, = atmospheric noise temperature 





2. RadioMet Operations 
UE Planner blocks: WFM 
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> Global scale forecast model: time and space boundary conditions to mesoscale models 
* RadioMetOP employs the ECMWF data: available every 12 hours, horizontal spatial 
resolution of 27 km and resampled at 18 km. 
» (Regional) Mesoscale forecast model: higher spatial resolution forecasts including explicit 
microphysics (global-scale model as an input) 
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Weather Forecast Module 


e RadioMetOP employs Mesoscale model 5 (MM®5) with runs set up to provide 1-to-5 day 
(24 to 120 h) forecast starting at 00:00 UTC for the 2012 
e ECMWF analyses provide the initialization of MM5 model and every next time step (6 
hours) ECMWF forecasts are used 
> WFM output: x,,,, on Cebreros area sampled every 1 hour 








2. RadioMet Operations 
UE Planner blocks: RPM 
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time series " 
ground-station specs 
Radio Propagation Module 
»  Radiative transfer models (RTM): RadioMetOP employs RTM based on modules of the 


Satellite Data Simulator Unit (SDSU) and Sky Noise Eddington Model (SNEM): 
* taking into account X tm, compute Loim and T tm, expressed as time series for several 


elevation angles. For the link budget, PDF and CDF curves of L,,,, and Tom are 
computed from the time-series 
> RPM output: PDF and CDF curves of Latm and Tatm computed from the time-series 








2. RadioMet Operations 
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Downlink Budget Module 


Data volume computation techniques (pass/sub-pass analysis): 
e Statistical: guaranteed weather availability (CD90%) 


e  Maximization: maximization of the average data return with a threshold on the lost 
frames (FIx,,=5%) 


DBM output: bit rate and minimum elevation angle per pass or sub-pass. 
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www 3. Data volume computation 


E 
Basic link budget formula xy O0 m {Crx A04), R, r, L (X am O) Toim (X aim (t)) 
0 


E 
Frame Error Rate computation FER(t) — frrr PL R, (0), X sy or 
0 
Transmitted frames per pass or sub-pass Fy (At) =— At 
lg 
R, 
Lost frames per pass or sub-pass Fy (At) = |, FER(t)at 
B 
Received frames per pass or sub-pass Fu UM) — FALCON) — Fy (At) 


R, = bit rate (constant in Af) discretized (BC discrete set of bit rates), rg = frame block 
lenght, Af= pass or sub-pass duration, 6,,,7 minimum elevation angle for the data 
transmission in A7 (210^). 
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3. Data volume computation: 


pass and sub-pass 


The Earth's rotation determines pass periods that are when there is a line 
of sight between the satellite and the ground station in Cebreros 
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During a pass period while the satellite rotates around Mercury it spins 
around itself 
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3. Data volume computation: 


=A parameters to optimize and 


antenna elevation profile [deg] 


benchmark 


Example of one pass with visibility interruptions and sub-passes for the year 2021 
80 rT T T T T T T T 


PASS 


70 | | | Parameters to optimize R,, Opn a 
^w higher 6,,,, allows higher R, (due to 


60! S | the lower slant path attenuation) with 
/ N the penalty of having a shorter A7. 
50 | | i | 
PA / N | S \ |, BENCHMARK (upper bound): ideal link 
/ s | budget calculation with a priori 
30 | p SUBPASSES e—— cu | knowledge of the time series of 
ua | atmospheric attenuation/sky-noise 
20 s. | , and the possibility of adapting the bit 
sub-passes \ rate every minute (analog bit rate), 


* blockages (visibility interruptions) 


10 ! i L 
2.422 2.423 2.424 2.425 2.426 2.427 2.428 2.429 
minute of the year 2021 " 10° 


Questions: 
* How can we take into account atmospheric random variability? 
* How to estimate the optimum bit rate and minimum value 6,,,,2 
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4. Rb and Om computation: 
statistical technique 


SAPIENZA 
Constraint: EJ/N? (E,/N,)) "H^ (E,/N;) tec? (E Nd mar 


BC mission: Coding scheme = Turbo Code ž > (E/N) r= 0.3 dB 
(EAN mar = 3 dB 
(E,/No) += 3.3 dB 


Statistical technique: R,,and 8, computation. Criterion: maximize Fz, 
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bod / Tatmp 


With «Fz4AT)» - Fr«( AT). 
Expected «F, ,(At)> and <Fp,(At)> are computed a posteriori after 9,. and &,.. determination. 
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4. Rb and Om computation: 
maximization Technique 


SAPIENZA 


Constraint: EJ/N? (E,/N,)) "H^ (E,/N;) tec? (E Nd mar 


BC mission: Coding scheme = Turbo Code ž > (E/N) r= 0.3 dB 
(EAN mar = 3 dB 
(E,/No) += 3.3 dB 


Maximization technique: &,, and @,,,,computation. Criterion: maximize «Foy? 


FER(t- fa [9 mj? Kyi L atm (x, (t). sf] — FER(At) = ET os Pr (5, ) dL, 
Oa 

] R 
[\e From BC < F yq (At) >= RI -—2.< FER(At) > | 
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5. Preliminary verification 
SAPIENZA e 
we chain (1/2) 


Preliminary veification chain to compute actual received and lost frames in a simulated 
operational scenario 


= ‘ = 
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Xatm: Laim: Vatm: Predicted weather conditions 
X am: Lain. V a! aC? ual weather conditions (for testing RMOP chain) 
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5. Preliminary verification 
SAPIENZA ° 
we chain (2/2) 


RMOP chain Verification chain 
Quantity Data source Quantity Data source 
Xatm (time series) MM5 X ary (time series) Meteo gauge 
Latms Tatm (time series) RIM L’ atm T atm (time series) SRM** 
RTM time-series/ 
Latm» Tatm Pdf, cdf) JPL tables* i j 


* long-term climatological statistics provided by NASA-JPL: monthly statistics 
obtained collecting radiometric measurements in Robledo for 19 years (1990-2009). 
** Simplified Radiative Model: Single rain slab with standard atmosphere conditioned 
to surface data 
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6. Results (1/3) 
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6. Results (2/3) 


DARTENZA 


Technique Database, timescale, analysis type y rary Mp 
JPL, Monthly, pass 964,9 100 

Statistical | RMOP, Daily, pass 985, 1 102,1 

(CD=90%) JPL, Monthly, s-pass 1146,2 118,8 


1189,6 1225.3 


JPL, Monthly, pass 
Maximization RMOP, Daily, pass 
(Frxij7590) JPL, Monthly, s-pass 


Benchmark - analog unlimited* Rb 2014,2 

Benchmark - analog limited** Rb 1728,8 

‘) Yearly received data volume “) percentage wrt reference case (in bold) 

* unlimited bit rate = no upper bound for R, 

** limited bit rate = upper bound for R, defined by the BC discrete set of bit rates 


Benchmark 





Monthly/daily timescale are referred to the timescale of the adopted statistics (cdf and 
pdf) in the statistical and maximization techniques. 
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6. Results (3/3) 
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Current: standard methodology for data return estimation, currently baseline for BC and 
using JPL data 

Advanced: natural evolution of the currently baseline for BC (still with JPL data) 
Potential (WF): methodology potentially applicable with weather forecast data 


. . a 
Statistical yearly rx DV yearly 





Technique & 
Methodology | | ob) mean 
Analysis type data Gbit % losses (% 
Current Statistical, pass Monthly-JPL 964,9 100 0,82 


Advanced Maximization,s-pass Monthly-JPL 1002,6 103,9 0.73 


Potential Maximization, s-pass Daily-RMOP 1205,5 124,9 0,76 
(WF) — 
(9 Yearly received data volume, (? percentage wrt reference case (in bold), © Weather 
Forecast : 


v 
Gain (in terms of percentage of received frames) between the WF -based 
approach and the considered methodology 
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7. Future works & 
acknowledgements 


SAPIENZA 


* Extend the period of analysis to more than 1 year (at least 5 years) 

* Reduce forecast time sampling from 1 hour to 10-15 min 

e Lack of experimental validation data: microwave radiometric 
measurements and radiosonde in Cebreros should be exploited 


* Extend the study to other frequencies (e.g. X-band) and to other 
ground station sites (e.g. Malargüe, Argentina) 
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